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Abstract 

This work studies heat transport of bond-disordered spin-1/2 chains. As an example, the XX case is analyzed, 
which corresponds to a model of noninteracting spinless fermions. Within the fermion representation, the single- 
particle eigenenergies are determined numerically, which allow one to compute transport coefficients. Since the 
ballistic transport properties of a clean chain are destroyed by disorder, the focus is on the frequency dependence of 
the thermal conductivity and on a qualitative comparison with the spin conductivity, both at finite temperatures. 
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Motivated by recent experiments on transport prop- 
erties of transition metal oxides [1] , intense theoretical 
work has recently been devoted to the study of heat 
conduction in one-dimensional spin-1/2 systems such 
as the Heisenberg chain, frustrated chains, and spin 
ladders [2]. Among these systems, the anisotropic 
spin-1/2 chain exhibits ballistic thermal transport 
properties at all temperatures T due to its integra- 
bility [3]. The thermal conductivity diverges in the 
homogeneous case, signaled by a finite thermal Drude 
weight [4]. Here, the usual decomposition of the real 
part of the thermal conductivity k into the Drude 
weight Dth and the regular part is utilized: Re k{uj) = 
DthS{iL!) + Kic^{uj), uj being the frequency. Ballistic 
transport is expected to be destroyed by randomness 
in the exchange couplings, resulting in a vanishing 
Drude weight. Therefore, one is interested in the reg- 
ular part Krog{u}), from which the dc-conductivity can 
be extracted by extrapolating to tj = 0. 

In this contribution, the effect of bond disorder 
on the thermal conductivity of XX chains is stud- 
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led numerically. This limiting case of the anisotropic 
spin-1/2 chain corresponds to free spinless fermions by 
means of a Jordan- Wigner transformation [5]. While 
particle transport in disordered fermion systems is a 
long-studied problem, closely related to the subject of 
localization (see, e.g., Ref. [6]), the thermal conduc- 
tivity K in these models has attracted less attention. 
Using bosonization, predictions were made for the 
concentration dependence of the thermal conductivity 
of spin-1/2 chains [7], and it could be very interesting 
to compare this to numerical results. Here, results for 
the frequency dependence of k are presented for the 
case of off-diagonal disorder at finite temperatures. 

The Hamiltonian in terms of spin-1/2 operators S"^'^ 
acting on site / and with periodic boundary conditions 
(PBC) reads: 



1 



(1) 



where N is the number of sites. It can equivalently be 
written in terms of fermionic operators cj : 



1 

H ^ Ji {cjci+i +h., 



(2) 
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In principle, a subtlety arises since the PBCs for the 
spin operators give rise to a nontrivial boundary term 
for the fermions, which depends on the number M of 

fcrmions [8]. This dependence, however, is neglected 
here, and PBC arc imposed for all M, since boundary 
effects can be expected to vanish for large system sizes. 

Introducing a spinor ip^ = (c|, . . . , cjy), the Hamil- 
tonian can be written as H = ip'^Aip, where ^ is a 
symmetric N x N band matrix, with nonzero elements 
= Ji/2 and Ai^n = Ji/2. Randomness in the J; 
is therefore called off-diagonal disorder, while a spar 
tially varying magnetic field realizes diagonal disorder. 

While a transformation to momentum space di- 
agonalizcs H for the translationally invariant case 
(J( = J), the computation of single-particle eigenen- 
ergies is still straightforward for random couplings 
(see, e.g., [8]). By means of a unitary transformation 
U with ci = Y2ij.^iiJ-Vfj. that diagonalizes the matrix A, 
the Hamiltonian can be written as H = ^nvjiVti i 
where are the single-particle eigenenergies. 

The energy current operator corresponding to 
Eq. (2) is J = j:i{clci+2-h.c.) = ij:^^ nlJ.'^n^- 

Within linear response theory, Kyee{ijj) is given by [6]: 

«reg(a>) = - E l^'-'^l' [/(^m) - /(e^)] ^('^ - Ae) , (3) 

where /(e) = l/(exp(/3e) + 1) denotes the Fermi- 
function, (3 = 1/T, and Ae = — e^. 

As an example, a Gaussian distribution of random 
couplings is considered: P{Ji) oc e"'"''""'' , with 
J — 1 and w — 0.2. The choice of the distribution does 
not affect the results, quantitatively consistent results 
for /t(w) are obtained with other Pi (e.g. box or binary 
distribution) by fixing the first moments of the dis- 
tribution function. Figure 1 shows results for Krcg(t<j) 
[pancl(b)] and, for comparison, arcg{(^') [pancl(a)], 
which is the spin conductivity corresponding to the 
current operator J^s = {i/2)'^i{c]ci+i — h.c). The 
parameters are T/J = 0.5 and = 1000, 5000. An 
imaginary broadening of 10^* was used. The overall 
form of both quantities is quite similar and resembles 
that known for cr reg(<^) [9]- Moreover, both curves 
exhibit a maximum at roughly the same frequency 
{uo/J « 0.015 for the parameters of Fig. 1). Both 
finite-size effects and statistical fluctuations are small 
in the low- frequency limit. 

The behavior of KTe%{ijf) at low frequencies is of par- 
ticular interest as it determines the dc-conductivity 
(see Refs. [10] for the case of clean spin systems). While 
the Drude weights vanish in the presence of disorder, 
the results of this work indicate finite dc-conductivities 
for both spin and thermal transport at finite tempera^ 
tures. This is illustrated in the inset of Fig. 1. A more 
detailed analysis of the finite-size scaling as well as the 
temperature dependence will be presented elsewhere. 
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(a) w=0.2,T/J=0.5 














|0.4 








) 0.0005 0.001 








1 1 

(b) 


\ 


1 — ^ — \ — ^ — 






- N=1000,#500 


' , 1 




— N=5000, #10 

1 



of , \ , \ , \ , 1 , 1 

0.02 0.04 0.06 0.08 0.1 



oVJ 

Fig. 1. Spin [panel (a)] and thermal conductivity [panel (b)] of 
an XX chain with a Gaussian distribution of random couplings 
J;, centered around J = 1 and with a width oi w = 0.2 (see 
text). Curves for N = 1000(5000) sites and 500(10) random 
realizations for the same parameters are displayed. An imag- 
inary broadening of 10~^ is used. The inset shows a blow-up 
of the low frequency region for the larger system size. 
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